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Abstract 


In this work, a design of transimpedance amplifier (TIA) for cryogenic scanning 
tunneling microscope (CryoSTM) is proposed. TIA with the tip-sample component 
in CryoSTM is called as CryoSTM-TIA. With transimpedance gain of 1 GQ, the 
bandwidth of the CryoSTM-TIA is larger than 200 kHz. The distinctive feature of 
the proposed CryoSTM-TIA is that its pre-amplifier is made of a single cryogenic high 
electron mobility transistor (HEMT), so the apparatus equivalent input noise current 
power spectral density at 100 kHz is lower than 6 (fA)?/Hz. In addition, “bias-cooling 
method” can be used to in-situ control the density of the frozen DX” centers in the 
HEMT doping area, changing its structure to reduce the device noises. With this 
apparatus, fast scanning tunneling spectra measurements with high-energy-resolution 
are capable to be performed. And, it is capable to measure scanning tunneling shot 
noise spectra (STSNS) at the atomic scale for various quantum systems, even if the 
shot noise is very low. It provides a powerful tool to investigate novel quantum states 
by measuring STSNS, such as detecting the existence of Majorana bound states in the 
topological quantum systems. 


1 Introduction 


High performance transimpedance amplifier (TIA) for cryogenic scanning tunneling mi- 
croscope (CryoSTM) is a key element [[ 2]. Recently, in Ref.[3], it was proposed a 
design of TIA for CryoSTM. TIA with the tip-sample component in CryoSTM is called 
as CryoSTM-TIA. In that design, the CryoSTM-TIA has a transimpedance gain of 1 
GO, a bandwidth of more than 200 kHz, and its equivalent input noise current power 
spectral density (PSD) that characterizes the circuit inherent noise is 21 (fA)?/Hz at 100 
kHz. The CryoSTM-TIA is capable of measuring the scanning tunneling shot noise spec- 
tra (STSNS) of various quantum systems at atomic resolution, even if the measured noise 
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current PSD is as small as a few (fA)?/Hz. In Ref.[3], it was illustrated how the CryoSTM- 
TIA is used to measure STSNS of quantum systems to investigate novel quantum states, 
such as detecting Majorana bound states (MBSs) in iron-based superconductors [4]. In 
the CryoSTM-TIA, cascade the pre-amplifier (Pre-Amp) and the post-amplifier (Post- 
Amp) to form an operational amplifier (OPA), and denoted as Macro-OPA. Pre-Amp in 
Macro-Amp is a differential amplifier made of a pair of high electron mobility transistors 
(HEMTs), where the HEMT is a kind of cryogenic GaAs MESFET [5} [6]. However, it is 
quite difficult to select cryogenic GaAs MESFETs with identical performances for pairing 
[2]. 

In this work, a design for the CryoSTM-TIA with a single HEMT in Pre-Amp is 
proposed. In the CryoSTM-TIA, cascade Pre-Amp and Post-Amp to form an inverting- 
amplifier (Inv-Amp), which is not OPA. The CryoSTM-TIA still has transimpedance gain 
of 1GQ and bandwidth more than 200 kHz. However, its equivalent input noise current 
PSD is only 5.3 (fA)?/Hz at 100 kHz, i.e. 1/4 of that in Ref.[B], since a noisy HEMT is 
reduced. Therefore, the lower tunneling shot noise of quantum systems can be measured 
at the atomic scale with this apparatus. Furthermore, only a single HEMT is used in 
Pre-Amp, avoiding the difficulty of matching two identical HEMTs. And, “bias-cooling 
method” can be used to in-situ reduce the inherent noise of the HEMT. With these 
advantages, this apparatus will be a powerful tool to investigate novel quantum states in 
various quantum systems, in the wider range of applications than that in Ref.[3], such as 
high-T, superconductors [7| [8], topological superconductors [4], MBSs [4] [9] [12], 
Andreev reflection [14], and Kondo effect [15], etc. 


2 Circuit of the proposed CryoSTM-TIA 


The circuit of the proposed CryoSTM-TIA is shown in Fig[1] It consists of several com- 
ponents: the single HEMT amplifier part of Pre-Amp shown in dashed box (al) of Figf1| 
the power supply of Pre-Amp in dashed box (a2), the post-amplifier (Post-Amp) shown 
in dashed box (b), the compensated feedback network shown in dashed box (c), and the 
signal source circuit shown in dashed box (d). The two stage amplifier made of Pre-Amp 
and Post-Amp is called as inverting-amplifier (Inv-Amp). Inv-Amp is connected with the 
feedback network to form TIA. The input G of TIA is connected to the signal source cir- 
cuit, to form the CryoSTM-TIA. The components placed in the cryogenic zone are shown 
in the dotted box. The parameters of all components of the CryoSTM-TIA circuit are 
listed in Table [I] 


2.1 Design of Pre-Amp 


As shown in Fig[] only a single cryogenic HEMT is used in Pre-Amp. The HEMT in 
this work is CNRS-HEMT (denoted as HEMT) [5] [6] developed by CNRS/LPN in France 
with the excellent cryogenic and noise characteristics. It is capable of operating under 0.5 
K with only 0.1 mW on the ideal operating point of “Vas = 100 mV and Ias = 1 mA”. Its 
parameters are listed in Table 1| where ex is its equivalent input noise voltage PSD and 
ie, is its equivalent input noise current PSD [6]. 

The single HEMT amplifier part of Pre-Amp is shown in dashed box (al) in Fig[1] 
The source of HEMT is grounded by a resistor R, and a small variable resistor Rgy;. 


© (al) 


Figure 1: Circuit of the proposed CryoSTM-TIA. Single HEMT amplifier part of Pre-Amp 
is shown in dashed box (al), power supply of Pre-Amp in dashed box (a2), Post-Amp in 
dashed box (b), compensated feedback network in dashed box (c), and signal source circuit 
in dashed box (d). The components placed in the cryogenic zone are shown in the dotted 
box. The parameters of all components of CryoSTM-TIA circuit are listed in Table [1] 


Rs = Rs + Rg, and a capacitor Cs of 0.1 mF is in parallel with Rg. Drain O of HEMT 
is connected to the load resistor R,. The resistor Rp is grounded through the variable 
resistor Rp1. Ry = Rp + Ry. The other end of Ry, denoted as Oz, is connected to the 
resistor Ra, and L as the other end of Ra is connected to P (the other end of Rı) through 
a small variable resistor Rp. Ri = Ro = Ry. P is connected to the power supply of 
Pre-Amp. P and L are grounded by the large capacitors C1 and Ca respectively. In the 
power supply, Rp = Ri + Ry, and Ry is a variable resistor. Ca is the input capacitance 
of Pre-Amp. HEMT, Rs, Rp, Ri, and R; are placed in the cryogenic zone, shown in the 
dotted box. Ci is the capacitance of the cable that connects gate G of HEMT and the tip 
in the CryoSTM. HEMT is as close to the tip as possible to reduce Cy, which can be less 
than 0.5 pF [18]. The operating point adjustment for Pre-Amp is shown in Sect [4] The 
Pre-Amp parameters are shown in Table [1] 

For Pre-Amp, when AC voltage V is applied to input G, the AC voltage difference 
between O, and O» is denoted as Vop- Ayp = Vop / V is the voltage gain of Pre-Amp. Since 
the gain-bandwidth-product of the CNRS-HEMT is gm /[27(Cgs+Cga)] ~ 1 GHz, the band- 
width of Pre-Amp is more than 30 MHz. In max{gm/(27Cs), 1/(27RsCs), 1/(27RgCo)} << 
f <3 MHz (ie. 30 Hz << f < 3 MHz with the parameters listed in Table [1), 


AvP © —9mfa, (2.1) 


where Ra = RL/(1 + Ruga). ie, Ayp ~ —20 in 30 Hz<< f < 3 MHz. The input 


Table 1: Parameters of all components of CryoSTM- 
TIA circuit 


HEMT (CNRS-HEMT) 


Gate-source resistance Ra >10 TQ 
Transconductance gm 40 mS 
Channel conductance ga 1 mS 
Gate-source capacitance Cs 5 pF 
Gate-drain capacitance Cga 1 pF 
Drain-source voltage Vas 100 mV 
Drain-source current las 1 mA 
er ((nV)?/Hz) 10 kHz 0.25 
100 kHz 0.07 
¿2 ((£A)?/Hz) 10 kHz 0.1 
100 kHz 1 
Pre-Amp 
Rs + Rsı 504502 Cs 11 yF 
Ry + Rui 150+5 Q Ri, R2 1kQ 
Rp 5+5 Q 
Cı 22 uF C2 11 uF 
PVR LM4050-10 R» 2 kQ 
Ri + Ra 4400 + 10 Q Vep +15 V 
Post-Amp with THS4021 as Rear-OPA 
Gao 97.5 dB fo 14.5 kHz 
Ca 1.5 pF Ra 1 MQ 
Rs 200 kQ Cii; Ci 150 pF 
Ra, Rea 100 Q Cri, Cre 50 pF 
VRp +15 V Ven -15 V 
Feedback network 
Rr 1 GQ CF ~ 3 pF 
Rk 300 kQ Ck ~ 0.2 pF 
Ce 10 nF 
Signal source circuit 
Ry >1 MQ Ci ~ 0.5 pF 
Note: + indicates the variable resistance range. Without 


specification, the default value after + is 0. 


capacitance of Pre-Amp is 
Ca = Cgs + (1 = Ayp)Cga- (2.2) 


ie., Ca = 26 pF in 30 Hz<< f < 3 MHz. The input resistance of Pre-Amp Ra is the 
gate-source resistance of HEMT, and Ra > 10 TQ, so it can be considered as infinity. 

The power supply of Pre-Amp is shown in the dashed box (a2) in Fig[1] A precision 
voltage reference (PVR) LM4050-10 is used in the power supply to provide a constant 
voltage of 10 Volt with typical temperature coefficient of +40 ppm/°C. The noises from 
PVR are eliminated by C1, Co, and Rr. 


2.2 Design of Post-Amp and composition of Inv-Amp 


In Fig[1] the Post-Amp circuit in the dashed box (b) is the same in structure as that in 
Ref.[3]. There is a commercial operational amplifier (OPA) in the circuit, called as Rear- 
OPA, which is the OPA with high gain-bandwidth-product, such as THS4021, OPA657, 


and LMH6624, etc. The Rear-OPA in this work is THS4021 [17]. The Post-Amp parame- 
ters are listed in Table[1] Ra and Ca are the input resistance and capacitance of Rear-OPA 
respectively. The feedback resistor Rf is connected to output O of Rear-OPA and its in- 
verting input. Rf is 200 kQ, and it placed in the cryogenic zone. Post-Amp also contains 
two cables that connect Q; and Q» in room temperture zone to outputs O4 and Oz of 
Pre-Amp in cryogenic zone respectively. The capacitance of the two cables is Ci, and Cig 
respectively, and they may vary from 50 pF to 150 pF. In this work, Cii = Cig = Ci = 150 
pF is assumed. 

A 1 meter cable with characteristic impedance of 50 2 has a distributed capacitance 
of 100 pF and a distributed inductance of 250 nH, and its distributed LC structure results 
in several resonant gain peaks in [50 MHz, 1 GHz] [20]. In Inv-Amp, Post-Amp and 
Pre-Amp are connected by the cables (Cable 0,Q, and Cable 02Q,), and the resonant 
gain peaks may cause the self-oscillations, as some unexpected electromagnetic couplings 
may exist among the outputs and the inputs of the amplifiers. I add two low-pass filters 
(Ra = Res = Re = 100 Q and Cri = Cra = Cr = 50 pF) with the upper cut-off frequency 
of 30 MHz at the Rear-OPA inputs, these resonant gain peaks can be greatly reduced, 
thus avoiding the self-excited oscillation. 

Cascade Pre-Amp and Post-Amp to form Inv-Amp. For the AC signal, the voltage 
gain of Inv-Amp is aa = aa(f). aa(f) can be expressed as 


aa = Ayp Ayr. (2.3) 


By the nodal analysis method, aa can be obtained. And then, with Ap expressed by 


Eq. (2.1), Avr can be obtained by Eq. (2.3). In max{gm/(27Cs), 1/(27RsCs), 1/(27RsC2)) << 
f <3 MHz (ie. 30 Hz << f < 3 MHz with the parameters listed in Table [1), 


Re 1 + 327 fRgLCir 1 
Rep 14+ 9207fRaCxr 14+ y janf Bre, : 


aa Rg, 


AvR x 


(2.4) 


where Ry, = Ry RL/(Ry +R), Cir = C¡+C;, and aa is the voltage gain of the Rear-OPA. 
In this work, Rp >> Ry and Ra >> Ru. And, a, can be approximately expressed as 
da = aa0/(1 + j f/f) in (0, 40 MHz] for THS4021. 

With the parameters in Table}1| the Inv-Amp performances are simulated by TINA- 
TI [19]. Figl2] show the curves of laa(f)lag and 4 (aa(f)) simulated by TINA-TI. By the 
nodal analysis method, aa(f) is calculated by the equations established with all compo- 
nents listed in Table |1| [21], which is also shown in Fig/2] The curves of laa(f)lag and 
Z (aa(f)) calculated by the nodal analysis method are identical with the simulated ones 


respectively. |aa(f)|ap and Z(aa(f)) are also calculated by Eqs.(2.1), (2.3), and (2.4), 


and the calculated curves in Fig[2]are almost identical with the simulated ones respectively 
in [3 kHz, 10 MHz], which verifies the correctness of Eqs.(2.1),(2.3), and (2.4). 


2.3 Frequency compensation of feedback loop 


In order to increase the bandwidth of the CryoSTM-TIA, for the high feedback resistor 
Rp with parasitic capacitance Cp, frequency compensation must be used in the feedback 
loop [1] [22]. In Figl] the compensated feedback network is shown in the dashed box (c). 
Rr is placed as near as possible to HEMT for reducing the capacitance between the gate 
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Figure 2: The Inv-Amp voltage gain aa(f). The solid curves are the curves of |aa(f)|ap, 
and the dashed curves are the curves of Z(aa(f)). The red curves are the TINA-TI 
simulation results, the green curves are the calculated results with the nodal analysis 


method, and the black curves are the calculated results by Eqs.(2.1), (2.3), and (2.4) in 
Article. 


G of HEMT and Ground. Taking Cù equalto kCp, where k is above 10°, adjust Ry equalto 
Rp/k, realizing RgCe = RpCp. The output voltage of TIA as V, generates the current Ip 
flowing to input G of TIA, so 


Vo a PR Rp 
~ jp 1+j2nfRkCk 14+ 520 f RC,’ 


where Ck is the parasitic capacitance of Rx 22]. Zp(f) can be considered as the 
impedance of the feedback network. In Ref.[22], it has been achieved in experiments to 
broaden the bandwidth of the feedback network with the very high feedback eee Rp 
of 10 GQ to MHz. In (0, 1 MHz], with the parameters listed in Table [1 E [Zr = 


Rp /|1 + j2r fRxCx| > Rp/1.008 and |Zp(f)| < Rp, so it can be considered that a 
equal to Rp. 
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Figure 3: At 4.2 K, Rp is 1.16 GQ and Cf is estimated as 3.36 pF. As Ri = 390 kQ and 
Ce = 10 nF, the black curve shows that the bandwidth of |Zp|/Rp is broadened to 905 
kHz, and the red curve shows zp as the phase of Zp is —31.5° at 905 kHz. 


In order to verify the above frequency compensation method, the following experiment 
is done. A resistor with 500 MQ(+25 ppm/°C) at room temperature is chosen as Rp, which 
is 1.16 GO at 4.2 K. Ry is a 380 kQ resistor in series with a0 ~ 20 kQ potentiometer, and Ce 


is a COG ceramic capacitor of 10 nF. In order to easily adjust the frequency compensation, 
a 3 pF capacitor is paralleled with the Rp, so that the Cp is about 3.2 ~ 3.4 pF (estimated 
as 3.34 pF). Fig.3 shows the experimental results with the frequency compensation, in 
which |Zp| = Rp/v2 at 905 kHz, and yzp as the phase of Zp is —31.5° at 905 kHz. 


2.4 Circuit stability of the proposed CryoSTM-TIA 


In Figh] the signal source circuit is shown in the dashed box (d), and its parameters are 
shown in Table The differential resistance of the tip-sample tunnel junction (TJ) is 
Ry, which is limited to no less than 107? Rp in this design. The capacitance of TJ is Cy, 
which is in parallel with Ry. And, Cy is estimated as several fF [3]. C = Ca + C1 + Cy in 
this work. Cy is at least two orders of magnitude less than Ca + Cr, so it can be ignored 
in C and C ~ Ca + Cr. The DC bias & modulated signal voltage source is denoted as 
BMS, which provides the DC bias V and sinusoidal modulated signal voltage V; for the 
CryoSTM-TIA. In the following simulation, C] is always taken as 0.5 pF. TIA connects the 
signal source circuit to form the CryoSTM-TIA. According to the parameters in Table [1] 
the performances of the CryoSTM-TIA can be simulated by TINA-TI. 
The loop gain Tj, of the proposed CryoSTM-TIA is 


Ti(f) = —aa(f)B(f) = aa (P/L, 
in which 6(f) is the feedback factor, and its reciprocal is, 


1/8(f) 1+ Zpll/Ry +1/Ra + ¡217 f(Ca + C1)]. (2.5) 


By Eq. (2.5), |1/8(f)| and < (1/6(f)) can be calculated. Fig/4{a) ) shows the calculated 
results of |1/G6(f)|ag and 4(1/B(f)) with Ry = +00, and Fig/4{b) ) shows those with 
Ry = 1 MQ. laalag and Z (aq) are also shown in Figl4] Both figures show |TLlag = 
laalas — [1/A(f)lan < —10 dB in f > 465 kHz and Z(T,) = Z(—aa) — (1/8) = 
Z (an) — 180° — Z (1/8) > —130° in f < 465 kHz. Therefore, the CryoSTM-TIA is stable 
with gain margin more than 10 dB and phase margin more than 50°. 


2.5 Voltage gain and transimpedance gain of the proposed CryoSTM- 
TIA 


With the compensated feedback network mentioned in Sect[2.3} it can be considered that 
Zp is equal to Rp in (0, 1 MHz]. Considering the TJ capacitance Cy, the TJ impedance 
should be Zz = Rj/(1 + j2afRyCy). And, C = Ca + Cr. As the AC input voltage V; is 
applied by BMS, the output voltage of the CryoSTM-TIA is Vo, and the voltage gain of 
the CryoSTM-TIA is A, = V,/V;. In (0, 1 MHz], by the nodal analysis method, Ay is 

Ae ie a (2.6) 


Z 1 Rp Rp Rp Cat Oy 
Jol aA aa Ry aa RA jenf 


Setting V; = 0 and applying a sinusoidal current source Å in parallel with TJ, the 
output voltage V, is generated at the output of the CryoSTM-TIA. A; = Vo Ji is called 
as the transimpedance gain of the CryoSTM-TIA. In (0, 1 MHz], A; is 

Rp 
Ay = (2.7) 
R R Rr(C uc 
1 = o TN j2rf El A D’ 
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Figure 4: TINA-TI simulation results for the voltage gain Inv-Amp aa, and 1/8(f) (a) with 
Ry = +00, (b) with Ry = 1 MQ. Both figures show |TLlag = |aalap — |1/8(f)lap < —10 
dB in f > 465 kHz and Z (TL) = 4(-—aa) — < (1/8) = Z (aa) — 180° — Z (1/8) > —130° 
in f < 465 kHz. Hence, the CryoSTM-TIA is stable enough. 


Disconnecting TIA with the signal source circuit, and applying a sinusoidal current 
source Ler into the input of TIA, the output voltage Vor is generated at the output of 
TIA. Ait = Vor/ Írr is called as the transimpedance gain of TIA. In (0, 1 MHz], Air is 


Lele ates — j2m f Sea! 


Considering Ra >> Rp and |aa| >> 1, Air is 


Rp 
A R C , 
1— Inf A 


Air = (2.8) 


In Fig[5} 5| |Air/Relap and Z (Air/Rp) calculated by Eq.(2.8) with the parameters in 
Table[1]is identical with their TINA-TI simulation results in E A 1 MHz], which verifies 
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Figure 5: TINA-TI simulation results for Air/Rp, Air is the transimpedance gain of TIA 
in the proposed CryoSTM-TIA. |Air(fnr)/Rrlap = —3 dB and Z (Air(fnr)/Rr) = 116° 
at for = 200 kHz. 


the correctness of Eq. (2.8) for Air. The upper cut-off frequency of Ajr(f)/Rp, i.e. its —3 
dB frequency is faer = 200 kHz. Comparing Eq. (2.7) with Eq. (2.8), as Ca + Ci + Cy & 
Ca +Cr ~ Ca and Rj > 107?Rp, the upper cut-off frequency of the CryoSTM-TIA is 
approximately equal to fpr. 


2.6 ‘Transient response of the proposed CryoSTM-TIA 


The simulation results for the transient response of the proposed CryoSTM-TIA are shown 
in Supplemental file 3 [24]. For the CryoSTM-TIA, the time taken from adding the input 
step signal voltage to the output response stably within a certain error is called transient 
response time tr. For the simulations, a resistor with a constant resistance Ro is instead 
of the tip-sample junction. When Ro > 30 MQ, for 0.1% error, tp < 5 us. When 5 
MQ < Ro < 30 MQ, for 1% error, tr < 5 us. When 1 MQ < Ro < 5 MQ, for 3% error, 
tr < 5 ps. 


3 Inherent noise of the proposed CryoSTM-TIA 


For the circuit of the proposed CryoSTM-TIA shown in Fig[1] the differential equivalent 
circuit with all noise sources is used to calculate its equivalent input noise. The details for 
the noise calculations are shown in Supplemental file 4 [28]. 


3.1 Equivalent input voltage noise and equivalent input current noise of 
Inv-Amp 


The equivalent input noise voltage and equivalent input noise current of HEMT are de- 
noted as ey and ¿y respectively. The resistors R¡, R2, and Ry are in the cryogenic zone of 
4.2 K, and their noises in f > 1 kHz are thermal noise, which can be neglected [3]. Ry is 
in the cryogenic zone of 4.2 K, the thermal noise voltage of the resistor Ry is ery. The 
thermal noise voltage of the resistor Re1 is e1, and that of Re2 is eg. The equivalent input 
noise voltage and equivalent input noise current of the Rear-OPA are denoted as e, and 
ia respectively. These noise sources are independent. The equivalent input noise voltage 
and equivalent input noise current of Inv-Amp are denoted as ea and 24 respectively. By 


the nodal analysis method and Wiener-Sinchin theorem, ignoring the minor terms, it is 
obtained that 


2 => 
E a A E (3.1) 
Ri, Alp | Ate | Ra Gen l 


= _ Ct R2 
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According to the parameters listed in Table |1| for Inv-Amp in this work, a= 0.28 


(nV)?/Hz and 7% = 0.14 (fA)?/Hz at f = 10 kHz, and eA = = 0.09 (nV)?/Hz and 7% = 4.9 
(fA)?/Hz at f = 100 kHz. For the Macro-OPA in Ref. [B], 1, eå = 0.5 (nV)?/Hz and a = 0.8 
(£A)?/Hz at f = 10 kHz, and e2 = 0.14 (nV)?/Hz and 1% = 20 (fA)?/Hz at f = 100 kHz. 
The resistor Ry in Inv-Amp in this work replaces the transistor H2 in the Macro-OPA in 
Ref.[3]. The noise generated by Ry is 2 orders of magnitude smaller than that generated 
by H2 in the Macro-OPA in Ref.[3], i.e. ee << Añez. Therefore, the equivalent input 
noises of Inv-Amp in this work are much lower than those of the Macro-OPA in Ref.[3]. 
In f > 10 kHz, e? = 2.25 (nV)?/Hz and 12 = 4 (pA)?/Hz [17]. In Eqs. (8-1), 


(3.2), and (3.3), exmRiL/Rú is one order of magnitude smaller than e2. In Eq.(3.1), 
(1 + RuL/Ra)*12/g?, is 2 orders of magnitude smaller than e?,. Further ignoring the mi- 


nor terms in Eqs.(3.1), (3.2), and (3.3), 


eh = ef + (ez + e + e2)/A?p, (3.4) 
= = i + + (27 f) Ch (e2 + ei + e2) /A?p, (3.5) 
eati = (iach) = —j2rfCa(el + e + €3)/A2p. (3.6) 


3.2 Equivalent input current noise of the proposed CryoSTM-TIA 

The equivalent input noise current PSD of the proposed CryoSTM-TIA is obtained as 
i =% + [L/R +1/RE + (2a f)?CR] el + 4kpT/ Re 

+ (1/Ry + j2nfCy) (eath + A/R) 


+ (1/Ry — ¿27 fC15) (Te + e/Rr) (3-7) 
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where Ciy = C1 + Cy [BB] [25] [28]. Putting Eqs. (3.4] (3-9, 6-5 , and (3.6) into Eq. (3.7), 
e, 0 + (1/Ry + 1/Re)* [ey + (+4 + a 


+ 4kBT/Rp + (27 f)? (Chek + C°(2 + ef + B)/Are| , (3.8) 
i.e. oe 
=i + (1/Ry +1/RpY A + 4kpT/Rp + i245, (3.9) 
where a= =(2177 (cas e+ 2CAC ¡(e + e? + e2) /A2p 


For the proposed co ee Rp = 1 GQ, Ca = 26 6 pF, CL = 0.5 pe, and Cy = 10 
fF. Rp and TJ are in the cryogenic zone at 4.2 K. As Ry = 1 MQ, i?, and its four 
components are listed in Table [2] The noise components of the CryoSTM-TIA in Ref. [3] 
are also shown in Table The CryoSTM-TIA proposed in this work has an equivalent 
input noise current PSD of 0.62 (fA)?/Hz at 10 kHz and 5.3 (fA)?/Hz at 100 kHz, which 


is much lower than that in Ref. [3]. 


Table 2: Noise components of CryoSTM-TIAs 


Pre-Amp type Single HEMT Amp Pre-Amp in Ref. [3] 
f (kHz) 10 100 10 100 
e? ((nV)*/Hz) 0.28 0.09 0.5 0.14 

Unit for the following terms is (fA)?/Hz 

a 0.14 4.9 0.8 20 
4kpT/ Rr 0.2 0.2 0.2 0.2 
Tas 0.0033 015 0.03 0.7 
(ata) 0.28 0.09 0.5 0.14 
i2 as total 0.62 5.3 1.5 21 


For the Epi-wafer of the CNRS-HEMT, as Si is doped in Al,Ga,_,As, each dopant 
gives rise to two types of electronic states: a shallow and delocalized donor level associ- 
ated with the normal substitutional site configuration, called as shallow donor state; a 
more localized acceptor level, known as DX” center, arising from a lattice distortion. In 
Al,Ga;_,As, as x > 0.22, DX” as a bound state is more stable than the shallow donor 
state. And, x = 0.37 for CNRS-HEMT in this work. There is a large repulsion barrier 
Ecap in the transition from the shallow donor to DX~, so DX” is frozen as T < 120 K 
[26 27]. 

According to the work on the noise mechanism of the CNRS-HEMT [5], “bias-cooling 
method” can be used to control the density of the frozen DX” centers in the doping area in 
the HEMT. This unique method, by which the in-situ modification of the HEMT structure 
is realized, can be used to reduce the noise generated by the HEMT. By cooling the HEMT 
with the positive gate-source voltage, more electrons in the HEMT are frozen into the DX” 
centers. At low temperatures, the absolute of the gate-source voltage required to maintain 
the HEMT at the ideal operating point is reduced, so that the gate leakage current is 
reduced, and thus the low frequency noises generated by the gate leakage current is reduced 
[5]. By the means, the inherent noise of the CryoSTM-TIA can be reduced. Moreover, 
if the HEMT used for the CryoSTM-TIA is not CNRS-HEMT, or its noise performances 
are not as good as those listed in Table [1] it can be remedied by “bias-cooling method”. 
This approach is very convenient for cases with only a single HEMT in Pre-Amp. 
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4 Operating state adjustment and DC tunneling current 
measurements 


For the CryoSTM-TIA in F ig[1] along with the circuit parameters listed in Table |1} its 
operating state can be adjusted as follows. (1) Disconnect Pre-Amp from Post-Amp and 
ground the input of Pre-Amp, i.e. gate G of HEMT. Keep Rp to 0. Adjust Rg, Rh1, 
and Rt, to achieve HEMT at the ideal operating point (Vas = 100 mV, Ig; = 1 mA) and 
the potentials at O} and Og equal. That is to say, Rg; and Ry, are adjusted to realize 
Ry = Rs + Vas/Ias = Rs +100 Q. (2) Cascade Pre-Amp and Post-Amp to form Inv-Amp, 
and input G is still grounded. Adjust Rp and Rt; to realize the DC output voltage Vom 
of Inv-Amp to 0 and HEMT still at the ideal operating point. (3) Connect output O and 
input G of Inv-Amp with the feedback network and disconnect input G from ground, to 
form TIA. Input G is still not connected with the signal source circuit, so the potential 
of input G and output O of Inv-Amp is still 0, since the input resistance of the transistor 
HEMT Ra can be considered as infinity. In the feedback network, Rp + Rk = Rp = 1 
GO. (4) Connect the signal source circuit to TIA, to form the CryoSTM-TIA. 

When f > 0, aa > aao, and aao is the DC voltage gain of Inv-Amp. |aao|ap ~ 
(gmRe/gaRu)ap ~ 94.5 dB, which is consistent with the simulation result laaolag = 92 
dB shown in Fig P| As the DC bias V; is applied by BMS, the DC resistance of TJ is 
R, and the potential at input G is Va, and the output voltage of the CryoSTM-TIA is 
Vo. Obviously, aaoVa = Vo, and (Vi — Va)/R = (Va — Vo)/Rr. The DC bias on TJ is 
V = V — Va, and the DC tunneling current I = (Vi — Va)/R. Is as an approximate value 
of I is 

ls = —Vo/ Ry. (4.1) 


And, the relative error is obtained as 
Er = |I; — I|/|I| =1/(1—aa0). 


As laaolag © 92 dB, Er < 30 ppm, consistent with the simulation results [29]. Since 
anoVa = Vo, the DC bias on TJ is 


"== (4.2) 


By Eq.(4.1), Vo = —Rels ~ —RpI = —RpV/R. By Eq.(4.2), Vi ~ V — RpV/(Raao). As 
the minimum of Ry is not less than 1 MQ, R > 1 MQ, so |Rp/(Rago)| < 1/40. Thus, 
V = V. With the measured V, I ~ J, and V ~ V; can be obtained. Therefore, the 
scanning tunneling current spectra J = I(V) (V €[VL, Vu]) can be obtained, where Vg 
and Vy are the lower limit and upper limit voltages for the measurements respectively. 

The drift of Von can be considered as the amplification of the drift of the Inv-Amp 
input offset voltage Vos, where Vom = aaoVos. The estimated drift of |Vos| is less than 
17 wV/°C [30]. The power of the power supply of Pre-Amp is less than 40 mW. The 
temperature control system based on the TEC devices can be used to control the 
temperature fluctuations of the power supply and those of Rear-OPA within 0.01 °C, so 
the fluctuations of Vos within 170 nV is guaranteed, i.e. the TIA output fluctuations 
within 170 nV. 

In order to prove the correctness of our proposed method for designing TIA, the N- 
Channel JFET (SST4393) is chosen to replace CNRS-HEMT, and TIA is designed and 
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fabricated with the same circuit design method and the same circuit topology. SST4393 
can work at 77 K. We measured the gain, noise, and transient response of JFET-based 
Inv-Amp and TIA at 77 K [32]. The experimental results are basically consistent with the 
theoretical calculations and simulations, which shows that the method of designing TIA 
in this work is correct. In addition, the designed JFET-based TIA can also be used for 
STM at 77 K, with larger bandwidth and lower inherent noise, comparing with traditional 
TIAs for CryoSTM (e.g. FEMTO DE-DLPCA-200 TIA). 


5 Applications for the proposed CryoSTM-TIA in spectra 
measurements 


For most applications, the modulus of the transimpedance gain of the CryoSTM-TIA 
|A;(f)| should be measured firstly. A;(f) is expressed as Eq.(2.7). Since laa (flag > 90 
dB in [1 kHz, 300 kHz] as shown in Fig(2| when Ry > 107°Rp, |Aj(f)| in [1 kHz, 300 kHz] 
can be approximately expressed as 


|Ai(f)| = Re/|1 — ¡27 fRp(Ca + Cr)/aa(f)]- (5.1) 


How to measure it has been described in Ref. [3]. 


5.1 Measurements of scanning tunneling differential conductance spec- 
tra by the proposed CryoSTM-TIA 


The differential conductance of TJ Gy = 1/Ry is the function of the voltage V applied 
to TJ. As the frequency f of the modulated signal voltage V; provided by BMS is low 
enough, such as f < 1 kHz, by Eq. (2.6), 
Rr 1 
A, = . e 
RyV) 1- OE 
Therefore, Rj(V) ~ [1/4, +1/aa(f)] Rp. When the measured |4,| < 1000, Ry(V) = 
Ry /|Ayl, since laa (flag > 85 dB in f < 300 kHz as shown in Figl2] With the measured 
A, | = |Vo|/|Vi|, the differential conductance spectra Gy(V) = 1/R3(V) (V €[VL, Vx)) can 
be obtained. 
Increasing the frequency f of the modulated signal can speed up the scanning tunneling 
differential conductance spectra measurements. In [1 kHz, 300 kHz], as Ry > 1 MQ, by 


Eq. (2.6), 


(AIS po * 
. [Z3(f)| |1 ¿27 fRp(Ca + Cr)/aa(f)| 

Therefore, 1/|Z3(f)| ~ |Av(f)|/|Ai(f)|, where |Aj(f)| is shown by Eq.(5.1). And, |A,(f)| 
and |Aj;(f)| can be measured. 1/|Z3(f)| = 1/R3 + (27 fCy)? can be obtained. Selecting 
two different frequencies fı and fa in [1 kHz, 300 kHz], |Z3(f1)| and |Z3(f2)| are obtained. 
Ry and Cy can be solved out from the measured |Z3(f1)| and |Z3(f2)|. 

Since the inherent noise of the CryoSTM-TIA is very small, the amplitude of the 
modulated signal voltage V; can be very small (< 10 uV), so that the energy resolution 
for STS is much improved. 
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5.2 Measurements of scanning tunneling shot noise spectra by the pro- 
posed CryoSTM-TIA 


The measurement method for the tunneling shot noise spectra by the proposed CryoSTM- 
TIA is basically the same as that introduced by Ref.[3], and it is only briefly described 
here. 

Before approaching the tip to the sample in the CryoSTM, Ry can be considered as 
infinity and Cy as 0. In this case, the output noise voltage PSD of the CryoSTM-TIA 
Ssu(f) can be measured, and the equivalent input noise current PSD i?(f) is 


P(f) = Seu A (5.2) 
where |A;(f)| is obtained by Eq. (5.1). 


To measure the scanning tunneling shot noise spectra (STSNS) of a quantum system by 
the CryoSTM-TIA, the distance between the tip and sample is adjusted, and the interval 
(VL, Vu] is selected, so that the shot noise measurements at V € Dy = {V| V% < V < 
Va, Gs(V) < 1 uS} are needed to study the physical properties of the quantum system 
[3]. In Dy, the output noise voltage PSD Ssum(f,V) can be measured. The tunneling 
noise current PSD is denoted as S}(f,V). The equivalent input noise current PSD of the 
CryoSTM-TIA is the function of Ry, so it is the function of V and denoted as i? (f, V). 
For Ssum(f, V), 


Sml V) = (V) BV| 141001. (5.3) 
With Eq. (5.2) and (5.3), it is obtained that 


sirv) = SL A) (5.4) 


From Eq. (3.9), 


AkgT eh 2 2,2 
Re + R? + (27 f)" | Cpe, + 2CAC1 


> => e+e +e 
RN=R+ E E 
vP 


2 “mn 


6(f,V) as the difference between i? (f, V) and i?(f) is smaller than 0.3 (fA)?/Hz in [10 
kHz, 100 kHz) and 0.1 (fA)?/Hz in[100 kHz, 300 kHz]. 9(f, V) can be neglected compared 
with Si(f, V), so long as the minimum of Gj(V) (V € Dy) is not too small. Therefore, 
Si(f, V) can be obtained as 


And i? ( f,V) is composed by the four parts as shown in Eq.(3.9). According to Table |2] 


Ssum V — Ssu 
sf, V) = ae w. 


by the measured Ssum(f,V), Ssu(f), and |Ai(f)|. And then, the SNTNS S]s(1) = 2Fel]1| 
can be extracted from S1(f, V). Here, J = I(V) (V € Dy), and I is the DC tunneling 
current as the bias V is applied. 

Compared with the CryoSTM-TIA in Ref.[3], the CryoSTM-TIA proposed in this work 
has the same bandwidth and transimpedance gain, but its inherent noise is much lower. 
The inherent noise of the CryoTM-TIA in this work is only 5.3 (fA)?/Hz at 100 kHz 


(5.5) 
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and its 9(f, V) < 0.1 (fA)?/Hz at 100 kHz, while that of the apparatus in Ref.[3] is 21 
(fA)?/Hz at 100 kHz and ô( f, V) < 0.28 (fA)?/Hz at 100 kHz. Therefore, as investigating 
novel quantum states of various quantum systems, the measurements performed with the 
apparatus proposed in this work is more accurate. 

For example, as investigating the existence of MBS in a magnetic flux vortex of the 
iron superconductor in CryoSTM [4], the tunnel junction resistance must be large enough, 
since the occurrence of incoherent Andreev reflection can be ruled out by weak tunnel 
coupling conditions [34]. In Ref.{4], the tunneling current I is quite low (see Fig.3(a) in 
Ref.{4] and Fig.s2(a) and (b) in its supplemental file). As the tunnel junction bias V is 
0.3 mV, I is only tens of pA, and the corresponding shot noise may only be only a few 
(fA)?/Hz. For the shot noise measurements in the system in Ref.[4], it is obvious that the 
CryoSTM-TIA proposed in this work is much more effective and accurate than that in 
Ref. [3]. 


6 Conclusion 


In this work, a design of transimpedance amplifier (TIA) for cryogenic scanning tunneling 
microscope (CryoSTM) is presented. TIA connected with the tip-sample component in 
CryoSTM is called as CryoSTM-TIA. The CryoSTM-TIA in this work has transimpedance 
gain of 1 GO and bandwidth more than 200 kHz, and its equivalent input noise current 
PSD is only 5.3 (fA)?/Hz at 100 kHz. In this CryoSTM-TIA, a single CNRS-HEMT 
Pre-Amp is instead of the differential Pre-Amp consisting of a pair of CNRS-HEMTs in 
Ref.[3]. The difficulty of matching the identical HEMTs is avoided. And, due to the 
reduction of a noisy HEMT, the apparatus inherent noise is only 1/4 of that in Ref.[3]. 
Furthermore, the inherent noise of the single HEMT in the circuit can be in-situ reduced 
by “bias-cooling method”. With this apparatus, the fast high-energy-resolution scanning 
tunneling spectra measurements can be performed and the very low tunneling shot noise 
of quantum systems can be measured at the atomic scale. This apparatus can be applied 
to investigate novel physical properties of various quantum systems, such as detecting the 
existence of Majorana bound states in the topological quantum systems. 
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Supplemental file 1: Voltage gain peaks for cables 


Figure s1-1 The circuit scheme for measuring the voltage gain peaks for cables. 


A 1 meter cable with characteristic impedance of 50 Q has a distributed capacitance of 100 
pF and a distributed inductance of 250 nH, which is shown as Fig. s1-1. In Fig. s1-1, R=1.5 kQ, 
the gain is Vo/Vin. The simulation results show its distributed LC structure results in several 
resonant gain peaks in [80 MHz, 1 GHz] shown as Fig. s1-2. We amplify Vo with Amplifier SA- 
220F5 (gain=200, fr=200 MHz), and then the gain peaks for the cable are measured with an 


oscilloscope. Fig. s1-3 show the measured gain peaks for cable. 


10.00 
0.00 
-10.00 
-20.00 


I (dB) 


-30.00 
-40.00 


-50.00 
0.00 


-1.00k 


HHY [deg] 


-2.00k 


-3.00k 
10.00MEG 100.00MEG 1.00G 


ME (Hz) 


Figure s1-2 The simulation results of the gain peaks for a 1 meter cable with a distributed 


capacitance of 100 pF and a distributed inductance of 250 nH. 
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Figure s1-3 The measured results of the gain peaks for a 1 meter cable with a distributed 


capacitance of 100 pF and a distributed inductance of 250 nH. 


Supplemental file 1: Voltage gain of Inv-Amp 


Figure s2-1 The Inv-Amp circuit with all components for Eqs.(s2.1-s2.8), and their values are 


shown in Table 1 in Article. 


With the nodal analysis method, the Inv-Amp voltage gain aa(f) is calculated by the 
following equations. All components for Eqs.(s2.1-s2.9) are shown in Fig.s2-1, and their values 


are shown in Table | in Article. 
(V,-V,)/R =(V.-V,)/R, +(e -V.)/Ra +P, A(R, + By) + JOC, (s2.1) 


joC,(V,-V,)=8,(¥,-%,)+8.(A-F)+(A-%)/R, +(V,-V,)/R, + joCy,, (s2.2) 


(P, 7.) /R, +(V, - Fp) /R, =V,/R,+ joC/,,, (s2.3) 
(V,-V,)/R, =(0, —V,)/R, + joc}, (s2.4) 
Em (VV) + ga (P V,)=V, [R + j0CY, + joC,,(V,-V,), (s2.5) 
V,=a,(V.-V), (s2.6) 


(V,-V.)/R, =+jaC,V., (s2.7) 


(Y, = V.)/R., = +j0C,V, > (s2.8) 


a, =V,/V,, (s2.9) 
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Figure s2-2 The Inv-Amp voltage gain aa(f). The solid curves are the curves of |aa(f)|as, and the 
dashed curves are the curves of Z(aa(f)). The red curves are the TINA-TI simulation results, the 


green curves are the calculated results with the nodal analysis method (i.e. Eqs.(s2.1-1.9)), and 


the black curves are the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. 


Fig. s2-2 shows the modulus of the Inv-Amp voltage gain laa(f)lag. The red solid 
curve is the TINA-TI simulation results, and the green solid curve is the calculated results 
with the nodal analysis method. They are identical with each other. The black solid curve 
is the calculated results by Eqs.(2.1), (2.3), and (2.4) in Article. In fe[3 kHz, 10 MHz], 


three curves are identical with one another. 


Fig. s2-2 also shows the argument of the Inv-Amp voltage gain Z(aa(f)). The red 
dashed curve is the TINA-TI simulation results, and the green dashed curve is the 
calculated results with the nodal analysis method. They are identical with each other. The 
black dashed curve is the calculated results by Eq.(2.1), (2.3), and (2.4) in Article. In fe [3 


kHz, 10 MHz], three curves are consistent with one another. 


Supplemental file 3: Transient response of the proposed CryoSTM-TIA 
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Figure s3-1 A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. TINA-TI simulation results for the transient response of the CryoSTM-TIA with the 
different Ro for the time interval of 10 us. The dashed curve is the step input signal Vi/Vins: for the 
different Ro. The solid curves are the output response V./ Vos. Here, Vost and Vins are their values 
respectively at 300 ms after applying the step input signal Vi, which are shown in Fig.s2-2. 


Transient response time ¢,<5 us. (b) is the zoom-in for the output response Vo/Vos in (a). 
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Figure s3-2 A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. TINA-TI simulation results for the transient response of the CryoSTM-TIA with 
different Ro for the time interval of 300 ms. (a) The step input signal Vi/Vins. (b) The output 
response Vo/Vost. Vo/Vost= 1 in t= 100 ms. (c) is the zoom-in figure of (b). 


A resistor with a constant resistance Ro is instead of the tip-sample junction for the 
simulation. Fig.s3-1 shows TINA-TI simulation results for the transient response of the 
CryoSTM-TIA with the different Ro for the time interval of 10 us. The dashed curve is 
the step input signal V;/Vins, and the solid curves are the output response Vo/Vos for the 
different Ro. Here, Vost and Vinst are their values respectively at 300 ms after applying the 
step input signal Vi. The error is definite as |Vo—Vost |/|Vost|. For the CryoSTM-TIA, the time 
taken from adding the input step signal voltage to the output response stably within a 
certain error is called transient response time f. As shown in Fig.s3-1 and Fig.s3- 
2(b)&(c), with Ro= 1 GQ and 30 MQ, for 0.1% error, <5 us. With Ro=5 MQ, for 1% 
error, 4< 5 us. And, With Ro= 1 MQ, for 2.5% error, tr<5 us. 


Supplemental file 4: Noise of the proposed CryoSTM-TIA 


S4.1 Noises of STM-TIA 
For Inv-Amp as shown in Fig.s1-1, the equivalent input noise voltage of the Inv-Amp is denoted 


as ea and its equivalent input noise current is ia, and their harmonic components of frequency f are 


E 
Ea and la respectively. iP can be obtained by the nodal analysis method. By Wiener-Khintchine 


A 


NA i EE, EI | 
theorem, cs can be obtained from A , [S4R1, S4R2]. The two matrix 
A A 


A A 


ine UN 
elements on the main diagonal are the equivalent input noise voltage PSD of the Inv-Amp e and 


its equivalent input noise current PSD if . The two matrix elements on the sub-diagonal are its 


equivalent input noise voltage-current PSD e,i, and equivalent input noise current-voltage PSD 


i A ‘ 

The Inv-Amp is connected to the feedback resistor Rr to form a TIA. The equivalent input noise 
voltage of the TIA is denoted as er and its equivalent input noise current is ir, and their harmonic 
components of fare Er and Jr respectively. The temperature of the feedback resistance Rr is T, and 


its noise voltage is ef and its harmonic component of f is Er. 


(b) 


Figure s4-1 (a) TIA circuit with the input short-circuit containing the equivalent input noise voltage of Inv-Amp 
Inv ea and its equivalent input noise current ia, and the output noise voltage of eoty ; (b) Noiseless TIA circuit with 
the equivalent input noise voltage of the TIA er as the input signal, and the output noise voltage of eotve ; the 
equivalency of the above two circuits means eoty = €otve - 
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For the TIA, the circuit containing all noise sources with the input short-circuit is shown as 
Fig.s4-1(a), and the output noise is eoty. The noiseless circuit with the equivalent input noise voltage 
of the TIA er as the input signal is shown as Fig.s4-1(b), and the output noise is éotve. For calculating 
the equivalent input noise voltage of the TIA, the equations are established on the equivalency of the 
above two circuits, 1.8. Coty = @otve. Therefore, by nodal analysis method, 


Er=EAa. 


Figure s4-2 (a) TIA circuit with the input open-circuit containing the equivalent input noise current of Inv-Amp 
Inv ea and its equivalent input noise current ia, and the output noise current of eori ; (b) Noiseless TIA circuit with 
the equivalent input noise current of the TIA ir as the input signal, and the output noise voltage of eotie ; the 
equivalency of the above two circuits means eori= €oTie - 

For the TIA, the circuit containing all noise sources with the input open-circuit is shown as Fig.s4- 
2(a), and the output noise is is eoti. The noiseless circuit with the equivalent input noise current of the 
TIA ir as the input signal is shown as Fig.s4-2(b), and the output noise is eotie. For calculating the 
equivalent input noise current of the TIA, the equations are established on the equivalency of the 
above two circuits, 1.8. €oTi =€oTie. Therefore, by the nodal analysis method, 


Ir=I1+HEa+ Ep/Rep. 


Hence, 
E 0 E 
T j= (E,+£,)+| ^ 
l I/R; L 
a A a 
From EA ai , by Wiener-Khintchine theorem, a Ertr can be obtained. The two 
tër Ly ie ï 


y E ir 


matrix elements on the main diagonal are the equivalent input noise voltage PSD of the TIA e2 and 


its equivalent input noise current PSD i2. The two matrix elements on the sub-diagonal are its 
q p T g 
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entr and equivalent input noise current-voltage PSD 
(s4.1) 
(s4.2) 


(s4.3) 


equivalent input noise voltage-current PSD 
(s4.4) 


2 * 
lrer: 


(s4.5) 


DEED 
e=e 
2_2 2 | p2 
2 =i +4k,T/R,+e/R2 , 
erip = (ne; ) *= ei, ten /R, 
Here, the noise voltage PSD of Rr is e = 4k,TR, . 


noise current PSD of a STM-TIA iz can be obtained [S4R1-S4R3] by 


The TIA is connected with the signal source circuit to form a STM-TIA. The equivalent input 


(YR + (2af)’ C3 es +R, + 2990, Jeri, +(I/R, — 72mfC, ize 


Here, Cy=Crt+Cj. Putting Eq.(s4.1), (s4.2), and (s4.3) into Eq.(s4.4), A is 


E 


a E, 


p 


(a2) 


P =i+ 
in =i, +4kaT [Re +e, [Y Ri + Re +(2nf) C3) 
+(1/R, + j2mfC, (ends +e, /Re) + (/R, - j2nfCy)(ige, +e /Re) 


Ro 


HEMT 


S4.2 Noises of the proposed CryoSTM-TIA 
Vp b Ro | 
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Figure s4-3 Circuit of the proposed CryoSTM-TIA. Inv-Amp of the proposed CryoSTM-TIA is shown in dotted 


box. 
For the proposed CryoSTM-TIA shown as Fig.s4-3 (1.e. Fig.1 in Paper). For the Inv-Amp, its 


e, iz, eri, , and ie, can be calculated by the nodal analysis method and Wiener-Khintchine 


theorem. Putting the results into Eq.(s4.5), the equivalent input noise current PSD of the proposed 
CryoSTM-TIA A is obtained. 


$4.2.1 The equivalent input noise voltage and the equivalent input noise current of Inv-Amp 


(a) 


(b) 


Figure s4-5 (a) Inv-Amp equivalent differential circuit with the input short-circuit containing all noise sources, and the 
output noise current of evay; (b) Noiseless Inv-Amp circuit with the equivalent input noise voltage of the Inv-Amp ea as 
the input signal, and the output noise voltage of @oave; the equivalency of the upper two circuits means @oay = @oAve. In (a) 
and (b), the triangle OPA is the Rear-OPA. 

For the Inv-Amp, the equivalent input noise voltage and equivalent input noise current of HEMT 
are denoted as en and iy respectively, and their harmonic components of fis are denoted as Ey and Jy 
respectively. gm is the transconductance of HEMT and ga is channel conductance of CH. The noise 
voltage of the resistors Ru is denoted as ern, and the harmonic components of fare denoted as Eru. 
The noise voltage of the resistors Ri, R2, and the feedback resistance Rr are denoted as ex1, er2 and er 
respectively, and the harmonic components of fare denoted as E 1, E12, and Er respectively. The noise 
voltage of the resistors Re and Re are denoted as e1, e2 and er respectively, and the harmonic 
components of fare denoted as E, E2, and Es respectively. Here, Ri= R2= RL. All noises of Ru, Ri, R2, 
Rei, Re2, and Re above 1 kHz are thermal noise. Ra = Ri//(1/ga) and Rut = Ru//Ri are denoted. The 


equivalent input noise voltage and equivalent input noise current of Rear-OPA are denoted as ea and 
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ia respectively, and the harmonic components of f are denoted as E, and /, respectively. These noise 
sources are independent of each other. For the Inv-Amp consistent of the Pre-Amp and Post-Amp, its 
equivalent input noise voltage and equivalent input noise current as ea and ia, the corresponding 
harmonic components of fare Ea and J, respectively. 

The equivalent differential circuit of the Inv-Amp containing all noise sources with the input 
short-circuit is shown as Fig.s4-5(a), and the output noise is @oay. The noiseless circuit with the 
equivalent input noise voltage of the Inv-Amp ea as the input signal is shown as Fig.s4-5(b), and the 
output noise is €oAve. For calculating the equivalent input noise voltage of the Inv-Amp, the equations 
are established on eoav = @oAve. Therefore, by the nodal analysis method, 

Ea ~ Ep + ErnZuv/(Avei Ru) ~Evi/(gmR) + Za(E2/Ri + El RO/Avpr 

+ Fa/Aypj + (E1-E2)/Avei — (1/2m + Zu/Ave;)La , (s4.6) 
where Aypi = gm/[1/Ritgat/(2nf/)Ci], Ave) = 2m/[1/Rit+getj2nf\Cr], Zar = Ruv/[1+ j(27nf)RurCir], and 
Cir =C; +C. 


IK 


Figure s4-6 (a) Inv-Amp equivalent differential circuit with the input open-circuit containing all noise sources, and 
the output noise voltage of eva; (b) Noiseless Inv-Amp circuit with the equivalent input noise current of the Macro- 
OPA ia as the input signal, and the output noise voltage of eoaie; the equivalency of the upper two circuits means 


€oAi= Codie. In (a) and (b), the triangle OPA is the Rear-OPA. 
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The equivalent differential circuit of the Inv-Amp containing all noise sources with the input 
open-circuit is shown as Fig.s4-6(a), and the output noise is @oai. The noiseless circuit with the 
equivalent input noise current of the Inv-Amp ia as the input signal is shown as Fig.s4-6(b), and the 
output noise is eoaie. For calculating the equivalent input noise current of the Inv-Amp, the equations 
are established on eoai= odie. Therefore, by the nodal analysis method, 

Ta~Tq HORA Car raza Ar) — J(27f/)CsaEui/(gmRi) + ¡OrRACA Zu (EL /RL+EY/RO/A we 

+ ¡Qr) (E1-E2)/Avei— ¡ORACAj EA ve ¡QC sam + CajZmL/A vrj)La, (s4.7) 

where Cai = Ces +Cgd 1+4wp1) , Csa = CestCga, and Caj = Ces tCgal 1+4 vw). 
For the measured frequency of fm < 300 kHz, 27fC;<0.3 mS with C¡=150 pF and 2xfC,<0.1 mS with 
C,=50 pF while 1/Ra= 1/Ri+ga=2 mS and 1/RuL = 1/RL+1/Ry=6 mS, so Avpi= 2m/[1/Ritgatj (2 nf)Ci] 
and Aypj= gm/[1/Ri+getj(2nf)Cir] are approximately simplified to -Ayp = gm/(1/Ri+ga) = gmRa, ZuL to 
Rut, and Caj to Ca. Thus, 


A (ate Belen 
= Jata UE y | e SEA 
A 1 0 i (2nf')C ApR i(2af)C, EaR j(2af)C, ) A RL (s4.8) 
1 Ra Ej 1 E-E, ahy 1+ Ru /Ra L 
JRf)C,) ApRe MICA) Ap I) (Ca + CR /Ra)) Em 


Ignoring the small quantities, such as the thermal noise of Ri, R2, and Rs, 


ne Ra )ean + É +e? te; é)|/4n+ (1+ Ra /R,) Re > (s4.9) 
oR ek 2 m Cte? +e? > a 

+(2 e+ (2 Ci 2122 4(2 C Ra c 2 4.10 

| wf) C ya ya SI 4 ( nf) Cx ya +(21f) at R ES a (s4.10) 


= Y. Re, tete R R \ P 
ei =h) = sevje E r` a 21+) C+ T a cae a > . (s4.11) 
d d m 


vP 


Here, e is the equivalent input noise voltage PSD of the CNRS-HEMT CH and 7 is their 
equivalent input noise current PSD. e? is the equivalent input noise voltage PSD of the Rear-OPA 


and P is its equivalent input noise current PSD. e? and e? are the thermal noise voltage PSD of 
Rea and Reoz respectively. 

e? =e2=4k,T.R, = 1.66 (nV)'/Hz, where Tr is 300K. e? =2.25 (nV)/Hz and i? =4 (pA)'/Hz 
in f> 10 kHz [S4R4]. In Eq.(s4.9)-(s4.11), Ri ff Ri ey is two orders of magnitude smaller than 


etete. (1+ Rap RY a e is one order of magnitude smaller than e in Eq.(s4.9). 
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(Cy, +C, Ra IR Y i g} is one order of magnitude smaller than Ci (e+e +e] / A in 
Eq.(s4.10). (Ca +C, Ru [ROH Ry, /R, 2/82 is one order of magnitude smaller than 


Cy (etete) / Az in Eq.(s4.11). Further ignoring the small quantities in Eq.(s4.9)-(s4.11), 


ee, t(ete tes) A, (s4.12) 
Y C (itita) A, (s4.13) 
ai hua) =-JQ2af)C, (2 te +4 |/4 (s4.14) 


$4.2.2 The equivalent input noise current of the proposed CryoSTM-TIA 
Putting Eq.(s4.12), (s4.13), and (s4.14) into Eq.(s4.5), the equivalent input noise current PSD of 


the proposed CryoSTM-TIA i is obtained as 


ER UR ARY EEEE) E er arca). 
(s4.15) 
where Cy = Cr+Cy and C= Ca+CI+C5. Eq.(s4.5) is Eq.(3.8) in Article. 


[S4R1] A. van der Ziel, Noise in Solid State Devices and Circuits, Wiley-Inter-Science, New York, 
(1986). 

[S4R2] Y.X. Liang, Low-noise large-bandwidth transimpedance amplifier for measuring scanning 
tunneling shot noise spectra in cryogenic STM and its applications, Ultramicroscopy, 234 (2022) 
13466. 

[S4R3] Z.H. Qian, Study on noise models algorithms and matrix descriptions for integrated circuits, 
J. Northeast Norm. Univ. 35 (2003) 41. 

[S4R4] Data sheet of THS4021 OPA, https://www.ti.com/lit/ds/symlink/ths4021.pdf.3. 
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Supplemental file 5: Estimating the DC tunneling current error for CryoSTM-TIA 
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Figure s5-1 For the proposed CryoSTM-TIA, Er=|f:-1//(1] vs. R (Re[1 MQ, 1 GQ]) is simulated by TINA-TI 
with V;= 1 mV (Black curve) and Vi=5 mV (Red curve). 


For the proposed CryoSTM-TIA, Fig. s5-1 shows the TINA-TI simulation results Er =|Js-/|/|J| vs. R 
(Re[1 MQ, 1 GQ]) with V¡=1 mV (Black curve) and V,=5 mV (Red curve). The simulation results 


show Er<30 ppm, consistent with the calculated results in Article. 
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Supplemental file 6: Estimation of Inv-Amp input offset voltage 


Figure s6-1 DC circuit of the inverting amplifier. 


Table s6-1 Parameters of THS4021 as the Rear-OPA in Post-Amp 


Open loop DC voltage gain aao 97.5 dB 
Input offset voltage drift |vord/°C 15 pV/°C 
Input bias current Jep, Zon 3 WA 
Input bias current drift |ipp|/°C, |Zon|/°C Estimated as 100 nA/°C 
CMRR 95 dB 
Ac = A20/CMRR 2.5 dB 


The Pre-Amp in the left dashed box of Fig.s6-1 is made of CNRS-HEMT CH. The source of CH is 
grounded via a variable resistor Rs. The source of CH is connected to Ri. Ri=R2=1 KQ, Rui =R1, Ri2= 
R2+Rp, and variable resistor Rp can be adjustable from 0 to 10 Q. Rui =Ri2=RL=1 KQ. Its input (the 
gates of CH) is grounded, and its two outputs O1 and O» are suspended. 

Adjust Rp, Rs and Ru, to let CH at the ideal operating point (Zas, Vas) =(1 mA, 100 mV), and the DC 
voltages at O is 0. The gate-source voltage Vgs = —JasRs, and Ve=Vas-Vzs. I assume Vgs =-50 mV for the 
ideal operating point, so Rs=50 Q and Ru= 150 Q. 

For the inverting amplifier, assuming a output voltage drift of vo for 1 °C, the voltage drift at the 
output O, is vı, and it at Oy is v2. The voltage drift at the output of the constant-current source is vs, the 
voltage drift at PVR is vp, and the voltage drift at the source of CH is vs. By the following equations, 
the input offset voltage drift of the inverting amplifier can be obtained. 


(vovi) RLH vv /RL= (vo v)/ Re, (s6.1) 
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(vvi) RL = ga(V1-Vs)-2mVs + ivp , (s6.2) 


(Vo-v2)/Ret (veva)/Ri= v2/Ru +ibn , (s6.3) 
vs/Rs = ga(Vi-Vs) —ZmVs 5 (s6.4) 
Aa0(V1tVott-V2) + Ac(VitVott +V2)/2 = Vo , (s6.5) 
Vos = Vo/@ao . (s6.6) 
aro ==gmR Y gaRh) . (s6.7) 


where, vos is the input offset voltage drift of InvAmp for 1 °C. From the above equations, vos and vs 
can be solved out as vos = ipp(1+gmRs)/gm7lngaRu/2m-Vp(1+emRs)/(2gmRt)-Vor( 1+2g4R1+gmRs)QgmRAL) 
and vs = vpRs/(2Rt)~ ibpRs. 

|vp[|<0.4 mV is the current drift of the constant-current source for 1 °C. Hence, |vs|<10 uV, and 
|vos|<17 uV. i.e., in Article, the drift of the gate-source voltage of CH in the inverting amplifier is less 
than 10 uV/*C, and the input offset voltage drift of the inverting amplifier is less than 17 uV/°C. The 
power of PVR is less than 40 mW, and the power of THS4021 is less than 200 mW. It is easy to control 
the temperature fluctuations of PVR and that of THS4021 within 0.1 °C with the temperature Control 


System, so the total input voltage drift can be controlled within 2 uV. 
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